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Abstract 

An off-line version of the Wall Interference Correc- 
tion System (WICS) has been implemented for the 
NASA Langley National Transonic Facility. The cor- 
rection capability is currently restricted to corrections 
for solid wall interference in the model pitch plane for 
Mach numbers less than 0.45 due to a limitation in 
tunnel calibration data. A study to assess output sen- 
sitivity to measurement uncertainty was conducted to 
determine standard operational procedures and guide- 
lines to ensure data quality during the testing process. 
Changes to the current facility setup and design rec- 
ommendations for installing the WICS code into a. new 
facility are reported. 


Nomenclature 


C D 

Non dimensional drag coefficient, 

Cl 

Nondimensional lift coefficient 

c x 

An output variable 

E 

Error 

M 

Mach number 

N 

Number of ports 

Q 

Dynamic Pressure 

u 

Uncertainty 

Ue 

Axial velocity determined by calibration 

u oo 

Free-stream velocity 

Ui,Vi,Wj 

Velocity perturbation; A r , \\ Z components 

X,Y,Z 

Cartesian tunnel axes 

a 

Angle of attack [rad] 
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8 Angle of sideslip [rad] 

A Correction 

e Blockage Factor 

a Population Standard Deviation 

a Population Standard Deviation Estimator 

Velocity potential 

Subscripts 

(%) Percent of value 

cow Corrected 

F Free air 

i incremental correction 

P Pressure 

T In tunnel 

t Total 

uric Uncorrected 

oc Freestream 

Superscripts 

(L) Confidence level 

Introduction 

The subject of wall interference has broadened im- 
mensely over the past century. Classical correction 
methods for this effect incorporate a mathematical 
boundary condition for the wind tunnel walls. Non- 
dimensional parameters have been tabulated and a 
number of formulas to correct spatially-varying wall- 
interference for aircraft-like models have been docu- 
mented in texts such as AGARDograph 109. 1 More 
often facilities are testing in regions that exceed the 
design specifications of the facility. Increased model 
size, for instance, introduces larger amounts of block- 
age and lift interference than the classical meth- 
ods can properly handle. Modern, more capable 
interference-correction methods, including boundary 
pressure methods were recently published in AGARD- 
ograph 336. 2 The advantage of the boundary pressure 
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methods over classical methods is their ability to re- 
spond to the actual conditions of the flow in the tunnel 
test section; additionally, boundary pressure methods 
are more complex than classical methods but can be 
implemented for online post-point/ post-scan applica- 
tions. 

Customers are also placing more stringent accuracy 
requirements on the acquired data. To reduce the 
uncertainties introduced in tunnel data due to the 
presence of the walls, thereby enhancing data quality, 
NASA Langley Research Center has implemented the 
Wall Interference Correction System (WICS) for solid 
wall testing in the National Transonic Facility (NTF). 3 
The code was originally developed by Ulbrich et al 4 9 
at NASA Ames Research Center for the 12-Foot Pres- 
sure Wind Tunnel. The WICS code uses a modified 
Hackett (Wall Signature) Method to determine wall 
interference corrections. This method was selected be- 
cause it is rapid and robust, and it has minimal impact 
on the facility in terms of instrumentation when com- 
pared to other boundary pressure methods. A brief 
description of the theory behind the WICS method 
and the implementation at the NTF is presented. 

The purpose of this paper is to determine code sen- 
sitivities to measurement uncertainty and to establish 
uncertainty boundaries. Guidelines for implementa- 
tion and operation are provided based on observations 
of this study. Although there is error inherent in the 
theoretical model due to limiting assumptions, this 
paper does not address the validity of the predicted 
corrections. 

Correction Accuracy Requirements 

Before a sensitivity study can be conducted on the 
WICS code, it is necessary to understand the accuracy 
requirements placed on wall interference corrections in 
general. According to Steinle and Stanewsky, 10 the 
requirements for transonic cruise are dominated by 
a one count (C D = 0.0001) accuracy for drag for a 
transport- type aircraft. Since the W ICS code is not 
used exclusively on transport-type aircraft in cruise, a 
more general requirement is needed. Newman and Ev- 
erhart 2 report the accuracy requirements listed below 
in Table 1. 


Type of Test 
High Lift 

Incremental 

0.2% C L kC D 

Absolute 

0.4% C L kC D 

Transonic 

1/2 count Cd 

1 count Cd 


Table 1 Generalized Industry Correction Accu- 
racy Requirements 


The Wall Interference Correction 
System (WICS) 

The following presents a brief overview of Ulbrich s 
extension of the Hackett wall signature method, in- 
cluding theory and the current implementation of 


WICS in the NTF. A more detailed description is 
given by Ulbrich et al 4 ' 9 and Iyer et al 3 The 
method applies a measured pressure boundary condi- 
tion which is the tared difference between the model 
installed condition arid the empty tunnel. The test 
article is comprised of singularities: a source-sink com- 
bination is used to represent the fuselage: sources are 
used to represent the viscous separation in the wake; 
and line doublets, typically placed on a lifting surface 
quarter-chord, are used to simulate the effect of lift. 
Additionally, a powered simulation can he represented 
by proper placement of a sink. Once the strengths 
of the line doublets and propulsion simulation sinks 
are known, the signature of these singularities is sub- 
tracted from the tared wall signature. The remaining 
signature is used to determine the strengths of the solid 
and wake blockage singularities. Hackett et al. 11,12 
originally used a “local” least squares fitting procedure 
of the wall pressure data, which was labor intensive 
because each pressure measurement used in the fit of 
the wall signature had to be individually selected for 
each data point. In an effort to improve the process, 
Ulbrich 9 introduced a “global” least squares fitting 
procedure that uses all available pressure ports. The 
use of panel codes to precompute normalized solutions 
for use in matching the signature was also introduced. 
Ulbrich further modified the Hackett method by us- 
ing balance measured forces and moments to calculat e 
the strength of the doublet singularities. The method 
uses tared data to assess the wall induced effects on 
the test article. The tared data are obtained by sub- 
tracting empty tunnel wall pressure port calibrations 
from the model-installed pressure data to remove first, 
order effects of the empty tunnel boundary layer and 
buoyancy. This taring of the measured data assumes 
that the additional second order change in the bound- 
ary layer displacement thickness due to the model 
in the tunnel is negligible. Normalized perturbation 
velocities are calculated from the subsonic potential 
equation using the method of images, for each type 
of singularity in discrete locations in the tunnel. By 
superimposing these model singularities with their cor- 
responding calculated strengths, interference velocities 
are computed. 

The blockage interference factor e is defined by 

€ =u, a -ur (,) 

where u f /L> = dipj/dx and v?/ « Incidence 

and sideslip corrections are defined by 


Wi , 

U e ' 

V e WT ~~ W F 

~ Coo u e 

(2) 

V < , 

U e vt — Vf 

(3) 

C, ' 

~ Coo U e 


where Uoo/U € is defined by 1+6. It is important 
to note that the integrated effect of the empty tun- 
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nel boundary layer growth and buoyancy due to im- 
proper wall divergence are already included in U € . The 
“empty-tunnel” calibration provides the perturbation 
velocity denoted by us , which can be used to correct 
the measured rnodel-in-tunnel perturbation velocity, 
xi’-p , by 

xir = u' T - u 6 

U e U c [4) 

The ratio uj /U e is used to determine the singularity 
strengths. 

General Implementation Requirements 

To implement the WICS code, 3 several tasks must 
first be completed. The first task is to ensure the facil- 
ity has an adequate number of static pressure orifices 
on the wall of the test, section. (The meaning of ad- 
equate will be discussed later in this paper.) Next, a 
perturbation velocity database (PVD) must be gener- 
ated using a panel method or the method of images. 
This PVD is used in the matching of the tared wall 
signature to determine singularity strengths, and to 
generate the field of interference velocities so that cor- 
rections can be computed. 

Once the orifices have been installed, an “empty tun- 
nel” calibration must be performed. This calibration is 
a function of independent test parameters ( e.g ., total 
pressure and Mach number) and support system atti- 
tude (e.g., pitch and sideslip angles). For semispan 
models, the empty tunnel calibration is the tunnel ge- 
ometry minus the model. For fullspan models the 
model support system is included, and the sting or 
support system kinematics must also be considered. 
Semispan and fullspan model installation details are 
provided subsequently. 

Code Output Results 

Figures 1-3 present a sample of the WICS code 
output from the NTF for a semispan model. Figure 1 
gives a row-wise comparison of the tared signature 
with the WICS-fit wall signature for a single test point 
at M — 0.2, a unr = 24° for a semispan model. The 
symbols represent measured, tared wall perturbation 
velocities and the solid line represents the global least 
squares prediction of the code. (One measure of va- 
lidity is how well these two match.) Mean primary 
corrections are displayed in Figure 2. These correc- 
tions are shown for an entire pitch polar (or group of 
data points). Blockage corrections (s) are averaged in- 
terference velocities along the fuselage centerline, and 
upwash (incidence) corrections (a,-) are averaged along 
the wing three-quarter-chord. Corrections to Mach 
number (AM) and dynamic pressure (A<2) are de- 
rived from the blockage factor. Coefficient corrections 
for the entire polar are plotted in Figure 3. Note for 
the test point of a = 24° in Figure 1 that signifi- 
cant corrections of AM - 0.0025 and a , = 1.4 are 
determined as shown in Figure 2. Correspondingly 


large corrections for the coefficients are also observed 
(Figure 3). Since the WICS code calculates the inter- 
ference corrections for specified reference planes in a 
volume surrounding the model, contours of the correc- 
tion variation can be plotted. 

Description of the NASA LaRC 
National Transonic Facility 

The NTF 13 is a fan-driven, closed-circuit, 
continuous- flow, pressurized wind tunnel, which 
is capable of testing at cryogenic conditions. The 
test gas is dry air for elevated temperature operation 
and nitrogen for reduced temperature operation. 
The settling chamber contains four anti-turbulence 
screens. A 15:1 contraction ratio entrance cone leads 
into an 8.2 feet square cross sectional test section with 
six inch triangular corner fillets which extends 25 feet 
in length then opens into a high spied diffuser. The 
operational envelope of the NTF encompasses a large 
range of test conditions. The facility can sustain a 
continuous airspeed from 0.1 to 1.2 in Mach number. 
Total pressure capabilities of the facility range from 
15 to 130 psi. The tunnel can operate at temperatures 
ranging from 150°F down to -320°F. These large 
ranges of conditions allow Reynolds number testing 
from 3 to 120 million per foot. NTF has the capa- 
bility to independently vary Mach number, Reynolds 
number, and dynamic pressure. 

Both fullspan and semispan model mounting sys- 
tems are available in the NTF. Fullspan models are 
supported by an aft mounted sting. This sting is 
attached to a vertically oriented arc sector used to 
change model pitch attitude over a range from -11° 
to 19°. The roll mechanism has a range of -90° to 
180°. Sideslip angles are achieved by using combined 
roll and pitch angles. A sidewall mounting system is 
used for semispan models. The angle of attack range 
is ±180°. The center of rotation for semispan and full- 
span models is tunnel station 13 feet. 

Currently, the NTF has 459 operational wall pres- 
sure orifices of which the WICS code uses 360. Looking 
downstream, Figure 4 shows the cross sectional dia- 
gram of the rows of pressure ports; it also serves as 
a map) between the NTF row numbering scheme and 
that used for WICS. Figure 5 unwraps the tunnel walls 
to show the current wall orifice layout. The filled port 
symbols are an example subset of port selections used 
for WICS. 3 

Approach and Rationale 

The original approach to define the WICS code sen- 
sitivities to measurement uncertainty was to use actual 
test data, perturb it using a known and controlled 
random error and bias, and evaluate the output. Fol- 
lowing this path is equivalent to declaring an existing 
data set to have zero uncertainty. This method was 
rejected because it gave no easy method for assessing 
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the magnitude of instrumentation error inherent in the 
current test data. Assessment of this error standard 
deviation is part of the ongoing Data Quality Assur- 
ance Program at LaRC conducted by Hemsch et al. u 

An alternate method was devised to obtain an ide- 
alized set of data from actual test data. The code, as 
part of the solution process, generates a global least 
squares fit of the wall signature. Additional coding 
was added to WICS to enable the discretization of t his 
wall signature fit. This result is reinput and a new fit 
of the wall pressure data is computed, thereby creat- 
ing a data set with no random error. This new wall 
signature and the corresponding test data are used as 
a baseline for the sensitivity study. In terms of the 
WICS calculation, force, moment, and reference con- 
dition data are point measurements. Therefore, only 
the wall pressure data were deemed necessary to be 
fitted (smoothed) to remove noise. The idealized data 
set assumes that there are no second order boundary 
layer or velocity perturbation effects, such as model 
induced nonuniformities in the tunnel wall boundary 
layer. Four test articles were used in this study; how- 
ever, all figures shown refer to a low aspect ratio, high 
lift, semispan model. A more detailed description of 
this analysis is given by Walker 15 including results 
from the other three test articles. 

The code has four major categories of error sources: 
the wall pressure signature, force and moment data, 
tunnel reference condition data, and the WICS math- 
ematical model. Given this information, the total 
uncertainty in the code output can be written as: 

Ut = /( U\\ PS , Up/M ? U ref ? Uw ICS M odet ) - (5) 

Each of the four categories in Equation 5 can be 
expanded. First, the wall signature is composed of 
precision and bias errors due to the wall pressure in- 
strumentation, thus yielding: 

Uwps — f{PwPSjBwps)‘ (6) 

Second, force and moment data consist of calculated 
lift, drag, and side force; roll, pitch, and yaw moments; 
and measured angle of attack. 

Up/M — f(Uip,U OP, Usf ? U PM , U YM ? u HM ) ■ (7) 

Third, the tunnel reference conditions, which in the 
NTF are determined from the Mach number reference 
system, will be examined. 

U REF = u MACH (8) 

Last, the WICS model incorporates the singularity 
representation, setup information, and the global least 
squares fitting process. The inherent error in the 
WICS model is beyond the scope of this paper and 
the corresponding uncertainty term is assumed to be 
zero or negligible for this study. A discussion of the 
remaining uncertainties is the subject of the remainder 
of this paper. 


Wall Pressure Measurements 

The NTF wall signature used for WICS computation 
currently consists of a total of 360 wall pressure orifices 
distributed in rows on the four walls. In Equation 6 
it was stated that the uncertainty in the wall pressure? 
signature could be broken into a precision error and 
a bias error. Precision error is defined as the random 
normal variation of a measured parameter about some 
expected value, whereas bias error is defined as the 
difference between some measured expected value and 
the “true value 11 of a parameter. The following sec- 
tions will evaluate the effects of these errors. As will 
be discussed, the magnitude of precision error is used 
to establish proper radial placement of rows, while the 
magnitude of bias error establishes the orifice distri- 
bution within the rows. 

Random Noise (Precision Error) 

From elementary statistics, a confidence interval is 
defined as the distance from the measured mean within 
which the true mean should lie given a specified confi- 
dence level or 

Cl (L) =x± (9) 
VN 

where x is the measured mean, b the measured stan- 
dard deviation of N samples or measurements and Z 
the standard normal statistic at a confidence or signif- 
icance level of L. Make note that when the number of 
samples is less than 30, Z is no longer standard normal 
but follows the student t distribution. 16 For simplic- 
ity, the standard normal distribution is assumed to 
hold for all values of A r . A confidence interval can also 
be thought of as the measured mean with some error 
band at a given significance level, written as 

x ± E {L) (10) 


If (j P is defined as the standard deviation of the wall 
pressure measurements, an expression for the error 
band of the wall pressures may be written 


E 


(L) _ 


Z^ l) ap 

y/N 


(ID 


Choosing a 99.72% significance level or a standard 
normal statistic of approximately 3 for data quality 
assurance purposes, the expected error band can be 
predicted for a range of op and N. Note that choos- 
ing a 99.72% significance level equates to one erroneous 
measurement out of approximately 350; whereas a 95%; 
significance level is equivalent to one out of 20. 

An estimate of the value of ap for the NTF is 0.002 
psi. This value was inferred from output of the tunnel 
on-line data system for low Mach numbers (0.1 to 0.3}. 


Random Seeding 

If the confidence interval is tight around the ex- 
pected value calculated from measured data, the true 
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mean is confidently and accurately predicted. For this 
case, at a fixed a, the predicted mean should be in- 
sensitive to random seeding of the random number 
generator used to introduce a prescribed measurement 
error. Increasing a in Equation 9 with fixed Z iL) and 
Ay will, at some value of a, cause the confidence* in- 
terval to become so wide that the mean cannot be 
predicted with accuracy. At this point the predicted 
mean will be increasingly sensitive to the effects of 
random seeding. 

Based on results presented by Walker, 15 the wall 
signature is insensitive to error due to random seeding 
for a p < 0.05 psi; beyond this, the wall signature can- 
not be accurately predicted. The global least squares 
fitting routine is able to adequately recover the un- 
derlying character of the pressure signature, thereby 
yielding a robust prediction technique. For values of 
op > 0.05 psi, the seed yielding the largest deviated 
case was conservatively selected to minimize an intro- 
duction of additional error in later analysis. 


Equation 12 may be written as 


E { ( !'J = uE ( p L) + b (13) 


Since this portion of the study is only concerned with 
E { p\ it follows that h — 0, which implies that when 
there is no error in the wall pressure system there is 
no error in the output variables. Thus, this model will 
only predict output, error due to random noise 1 in the 
wall signature. Using Equation 11, Equation 13 may 
be re-expressed as 



= a 


Z^'^op 

y/N 


Solving Equation 14 for N yields 


(14) 


N = 


aZ iL) a P 
El L) 


(15) 


Variation of a P 

Using the typical value of a P as a guideline, six val- 
ues, ranging from 0.001 to 0.1 psi, encompassing three 
orders of magnitude were used to define the code sensi- 
tivity to variation of a P . The wall signature is not very 
sensitive to random perturbation as discussed in the 
previous section. The insensitivity of wall pressure to 
random perturbation translates, as would be expected, 
into output insensitivity as evidenced in the incremen- 
tal blockage results presented in Figure 6. This figure 
is representative of the sensitivities of each output vari- 
able examined over the data run for various values of 
op. The maximum and average values for the data run 
are shown. Figure 6a shows the scatter of incremental 
blockage error for the entire data run, and Figure 6b 
shows a linear-logarithmic plot of the maximum and 
mean values for each value of a P . The solid and dashed 
lines represent a linear fit of the maximum and average 
values for each value of a P , respectively. This figure is 
representative of all the output and on this logarith- 
mic plot it shows a strong positive linear correlation 
between the absolute and average error in the output 
variables and a P . It is also evident that the distance 
between the fit of the maximum and mean errors be- 
comes larger as (t p increases, demonstrating that the 
scatter in the error becomes large with relatively large 
o P . 

Variation of N (Loss of Pressure Data) 

A precision-sensitivity model that will enable the 
prediction of output error due to random noise in the 
wall signature can be developed as follows. It is un- 
derstood that 

= /(44 (i2) 

where C x is any output variable. If / is assumed linear 
due to the strong linear correlation as seen in Figure 6, 


Equation 14 allows ) to be predicted for variation 
in N. and Equation 15 allows for the prediction of the 
number of required pressure ports for a given a P , Z^ L \ 
and E ( r L) . 

By regressing the maximum values of error for each 
of the output variables, a conservative estimate of the 
value of a for Equation 13 may be determined for each 
output variable. Table 2 tabulates these values. Us- 
ing Table 2 and Equation 15, prediction plots at the 
99.5% confidence level were made corresponding to 
several of the output variables. Figures 7 and 8 show 
the predicted error for the primary corrections e and 
a n respectively. In these figures, the predicted num- 
ber of ports required is plotted against a P for lines 
of constant error in the output variables. This allows 
the determination of a required number of ports for a 
given ap and a defined maximum error level. Make 
note that the units of each error level correspond to 
the units of the output variable. As an example, for 
the NTF typical ap = 0.002 psi, the number of ports 
required to ensure incidence correction error predic- 
tion below 0.005° is 240. The corresponding number 
to ensure that the blockage correction is missed by no 
more than 0.01%. is 120. Accuracy of the angle of at- 
tack is obviously the dominant criteria for selecting 
the number of ports based on error in primary cor- 
rections, while blockage requirements can be achieved 
with a relatively small number. In accordance with 
the generalized requirements placed on testing in Ta- 
ble 1, Figure 9 gives the predicted error in terms of 
ACo t , and Figures 10 and 11 give the predicted error 
in terms of A C L{%) and A C Dt{%) , respectively. 

Row Placement 

The eight semispan and ten fullspan wall pressure 
orifice configurations defined in Tables 3 and 4, re- 
spectfully, were used to evaluate output sensitivity due 
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Variable 

a 

£ 

0.1887 

AM 

0.0382 

AQ (psf) 

129.51 

a i: 

2.G32 

AC l 

0.0118 

^Cl(%) 

1.G52 

AC m 

0.0167 

AC d 

0.0349 

AC d , 

0.0317 

AC n , 

0.0382 

AC Dt %) 

14.065 

D b (%) 

17.461 

ACd,(%) 

14.868 


Table 2 Regressed Values for the Coefficient a 
in Equation 15 for Prediction of Error in Output 
Variables due to Variation in a p and N 


Case 

Name 

WICS Rows 

N 

1 

FulLSig 

1-5,7-12 

260 

2 

Alternating 

1,3,5,7,9,11,12 

161 

3 

Mirror.CLine 

5,7,12 

67 

4 

TunneLCLine 

3,9,12 

64 

5 

Far -Field 

1,11,12 

62 

6 

Min.Top 

7 

26 

7 

Min_Bott 

5 

25 

8 

Min -Side 

12 

16 


Table 3 Semispan Row Hookup Configurations 


Case 

Name 

WICS Rows 

N 

1 

FulLSig 

1-12 

290 

2 

Alternating 

1,3,5-7,9,11,12 

191 

3 

FourCorner-Side 

1,5-7,11,12 

143 

4 

FourCorner 

1,5,7,11 

97 

5 

Tunnel.CLine 

3,6.9,12 

94 

6 

14x22-Sim 

6.9,12 

72 

7 

Min_Far_Side 

6 

30 

8 

Min.Top 

9 

26 

9 

Min_Bott 

3 

22 

10 

Min_Near_Side 

12 

16 


Table 4 Fullspan Row Hookup Configurations 


to row placement configurations and to verify predic- 
tions in Figures 7 - 11 at various values of op. These 
configurations include the full measurement, near field 
(radially), far field (radially), and single row resolu- 
tions of the wall signature. In addition, approxima- 
tions to other facility configurations were evaluated. 
Walker 15 showed that the error predicted in Figures 7 
- 11 is a fairly conservative estimate with notable ex- 
ceptions due to orifice placement. 

In the random perturbation of the wall signature, 
the doublet strengths which define the lift influence are 


not affected. Therefore, all randomness is transferred 
to the blockage component of the pressure signature- 
arid the blockage factor will be more sensitive to this 
type of experimental noise. The error in blockage in- 
creases with the decrease in the number of active ports, 
which corresponds to a less accurate fit of the singular- 
ity strengths. Error in the blockage factor was found 
to be approximately the same order of magnitude for 
numbers of ports greater than or approximately equal 
to 150. Row configurations with less than 150 ports 
are therefore less desirable. Since the upwasli correc- 
tion incorporates both blockage and lift interference, 
the overall trends in the upwash correction factor will 
differ slightly from those of the blockage. As with the 
blockage correction, the magnitude of error in the up- 
wash correction tends to be approximately the same 
for cases in which the number of ports is greater than 
approximately 150. This again suggests that row con- 
figurations with less than 150 ports are inadequate, 
unless special care is taken in the placement of rows 
and orifices. It is also important to note that the mea- 
sured pressure signature from rows of orifices in the 
spanwise far-field of the model yield a higher level of 
error than those rows positioned near the model. This 
is primarily due to the radial distance of the ports from 
the singularity representation of the model. Given 
these trends, it is recommended that the majority of 
the orifices be located in the vicinity of the model. It 
is noteworthy to point out that a higher average error 
exists for measurements from rows of pressure orifices 
placed directly below the model as opposed to directly 
above the model. This suggests a slight favoring of the 
ceiling orifices, and in terms of maintenance and orifice 
integrity, ceiling pressure ports are favorable. Rows of 
pressure orifices in the spanwise far-field and small to- 
tal numbers of orifices yield less accurate results due 
to the inability of these configurations and numbers 
of ports to adequately resolve the wall signature. This 
results in the inability of the method to accurately dif- 
ferentiate between the lift and blockage interference in 
these instances. While some of these rows are less im- 
portant for pitch-plane corrections, the sidewall rows 
will become essential when full three-dimensional cor- 
rections are implemented. 

Bias 

Bias error was introduced into the wall pressure sys- 
tem by adding a constant amount to each wall pressure 
measurement. The instrumentation used to record 
wall pressure data has a quoted accuracy of ±0.0025 
psi; additionally, bias of this magnitude was observed 
during the acquisition of data used for this study. This 
range is similar to the precision error case; therefor e, 
the same three orders of magnitude study was con- 
ducted to determine sensitivities of the code to bias 
error. Two separate studies related to bias error were 
conducted: a global biasing (i.e., the entire wall signa- 
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ture) and selective biasing of chosen pressure modules. 

Global Biasing 

Biasing of the wall signature on a global scale has 
the same effect on blockage as changing the Reynolds 
number; when bias is increased the measured blockage 
displays an apparent decrease. When Reynolds num- 
ber increases, the tunnel wall boundary layers thin and 
blockage decreases. This suggests that bias error in the 
wall pressures can manifest as a pseudo- Reynolds num- 
ber effect. Figure 12 shows the biased wall signature 
for two values of Bp as compared to the baseline case. 
This figure shows the biased wall signatures shifted 
from the baseline case and gives the corresponding 
fit of each case. As demonstrated by Walker, 15 the 
blockage is affected significantly by bias error; whereas 
the upwash correction contains only a slight variation. 
Figure 13 shows the sensitivity of the primary correc- 
tions to an increase in B P . In contrast to the random 
noise study, a minimum error now exists because the 
measured values have been shifted away from the true 
solution. Very little difference in the blockage is seen 
in Figure 13a, but the upwash correction in Figure 13b 
contains a significant difference between the maximum 
and minimum values of error for each level of bias. 
This difference in the upwash is due to the fact that the 
upwash correction contains the superposition of both 
lift and blockage interference, and any differences seen 
in the blockage factor are magnified in the upwash cor- 
rection. Figure 14 shows the error in terms of the force 
coefficients in percentage of the corresponding uncor- 
rected coefficient, and Figure 15 shows the change in 
the coefficients. These figures can be used to deter- 
mine the level of bias allowed according to the error 
levels in the general requirements for testing (see Ta- 
ble 1). Error in both lift and drag corrections has a 
strong linear correlation with Bp . 

The effect of the bias is clearly seen in the following 
example. Several repeat data runs were done dur- 
ing the test of a low-aspect ratio, high-lift, semispan 
model. Normally, before each set of data runs the wall 
pressure system was fully calibrated while the tunnel 
was offline. Calibrations performed during a run are 
zero calibrations. A full calibration entails fitting a 
fourth order curve through five data points dispersed 
over the testing range. During a zero calibration the 
curve is assumed to have shifted and only the constant 
or zero term is adjusted for offset, the other coefficients 
remaining constant. For one particular run the refer- 
ence pressure drifted when the tunnel came online and 
only a zero-calibration of wall pressure s} 7 stem was con- 
ducted prior to taking data (the full calibration of the 
wall pressures was performed while the tunnel was of- 
fline). For this case it resulted in a biasing of the wall 
pressure data. During the same data run another zero 
calibration was performed which unbiased the data. 
Corrections were computed for these data. Figure 16a 


shows the blockage results for the three repeat data 
runs. The step change in the data at a « 14° plotted 
with deltas are the results where the zero calibration 
was performed during the run. After calibration the 
offset blockage results collapse and follow the trends of 
the two other repeat runs. Figure 16b shows three re- 
peat runs where a full calibration was performed offline 
with no zero calibration when the tunnel came online 
and no problems were observed. Although this bias- 
ing phenomena can be duplicated with artificial error, 
further investigation should be done into the calibra- 
tion process of wall pressures. The zero calibration 
taken before acquiring data was probably done before 
the static reference pressure reached equilibrium and 
thereby biased the data. When the second zero cal- 
ibration was performed the static reference pressure 
had probably equilibrated and the bias was then ad- 
justed out. 

Another significant outcome of the global biasing of 
the wall pressure data is that of port placement for 
testing. By subtracting the fit of the wall pressures 
from the biased wall pressures, the ports which exert 
the most influence on the wall fit can easily be seen. 
This effect can also be seen in Figure 12 by comparing 
the biased data with its corresponding WICS fit. The 
areas where the fit has relatively large deviation from 
the data are upstream of the model. The data from 
the ports in these areas of large error had an insignifi- 
cant effect on the fit generated by the code. Therefore 
an optimal port configuration for using the WICS code 
would be to have a vast majority of the pressure ports 
in the vicinity of the model and approximately one 
fuselage length downstream where the wake is influen- 
tial. The primary reason for this is the distance from 
each port to singularity representation dictates how 
much influence that port has on the final output. This 
also suggests that tunnel aspect ratio will have some 
effect on the fitting of the wall signature since the dis- 
tance from the model representation singularities to 
the pressure ports may be larger or smaller. 

Selective Biasing 

Selective biasing was performed by adding a pres- 
sure increment to all the data from select pressure 
sensing modules. It can be inferred from results pre- 
sented by Walker, 15 that error due to biasing becomes 
smaller with fewer ports being biased so that the bias 
is in effect averaged out over the other ports. This is 
seen predominately in the error of the blockage factor, 
with one notable exception that the upwash correc- 
tion seems to be more sensitive to a biasing of the 
floor pressures than the ceiling. 

Force and Moment Measurements 

Of all the forces and moments measured on the 
model, lift is the primary input to the WICS code. 
As stated previously, the lift force is used to calcu- 
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late the strengths of the line doublets. This means 
that error in the lift force will cause more error in 
the upwash correction than blockage. Figure 17 ev- 
idences this statement by showing resultant error in 
the primary corrections due to introducing error in lift 
according to the generalized industry requirements in 
Table 1. Here, the blockage and upwash correction er- 
rors, Sc and Sctj are plotted against o U nc for error in 
measured lift of 0.2% and 0.4%. These lift errors are 
insignificant lor blockage; however, they can achieve as 
rmirh tls 50% of the allowable error for angle of attack 


U f /m and U ref terms become negligible and drop 
from Equation 5. This allows the total uncertainty to 
be approximated by Equation 6. Then 

U t ~ Uwps = f(Pwps,B\vps) (1G) 

The combination of errors is done according to the 
root sum square method (RSS). With the use of this 
method, any error less than half the size of a larger er- 
ror is negligible. According to results presented in the 
wtinn n n wall nres.su re measurements, the amount 


13 Fuller, D., Gloss, B., and Nvstrom, D., “Guide for Users of 
the National Transonic Facility,” NASA TM 83124, July 1981. 

l4 Hemsch, M., Grubb, J., Krieger, VV., and Cler, D., “Lang- 
ley Wind Tunnel Data Quality Assurance- Check Standard Re- 
sults,” AIAA Paper 2000-2201, 2000. 

15 Walker, E., Sensitivity of the Wall Interference Correc- 
tion System to Measurement Uncertainty , Master’s thesis, The 
George Washington University, 2000. 

16 Coleman, H. and Steele, W., Experimentation and Uncer- 
tainty Analysis for Engineers , Wiley-lnterscience, 1989. 
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Fig. 1 Wall Signature of a Model at a unc = 24° 
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(Looking downstream) 



Fig. 4 Cross-Sectional Pressure Measurement Row Setup of the NTF 
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Fig. 5 Wall Orifice Layout for NTF 
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(Open symbol =» Port is hooked up for this test Symbol + => port is bad or not hooked up Filled symbol => port selected for WICS) 
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Fig. 9 ACd ( Error Prediction 
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Fig. 12 Wall Signature Comparison for Bp — ±0A psi 
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Fig. 14 Percent Force Coefficient Correction Sensitivity to Variation in B P 


24 OF 27 


American Institute of Aeronautics and Astronautics Paper 2001-0159 



I5AC 
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Fig. 15 Force Coefficient Correction Sensitivity to Variation in D , 
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